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Although the biochemical pathways for the anaerobic degradation of many of the 
hydrocarbon constituents in petroleum reservoirs have been elucidated, the mechanisms 
for anaerobic activation of benzene, a very stable molecule, are not known. Previous 
studies have demonstrated that Geobacter metallireducens can anaerobically oxidize 
benzene to carbon dioxide with Fe(lll) as the sole electron acceptor and that phenol 
is an intermediate in benzene oxidation. In an attempt to identify enzymes that might 
be involved in the conversion of benzene to phenol, whole-genome gene transcript 
abundance was compared in cells metabolizing benzene and cells metabolizing phenol. 
Eleven genes had significantly higher transcript abundance in benzene-metabolizing cells. 
Five of these genes had annotations suggesting that they did not encode proteins that 
could be involved in benzene metabolism and were not further studied. Strains were 
constructed in which one of the remaining six genes was deleted. The strain in which 
the monocistronic gene Gmet 0232 was deleted metabolized phenol, but not benzene. 
Transcript abundance of the adjacent monocistronic gene, Gmet 0231, predicted to 
encode a zinc-containing oxidoreductase, was elevated in cells metabolizing benzene, 
although not at a statistically significant level. However, deleting Gmet 0231 also yielded 
a strain that could metabolize phenol, but not benzene. Although homologs of Gmet 
0231 and Gmet 0232 are found in microorganisms not known to anaerobically metabolize 
benzene, the adjacent localization of these genes is unique to G. metallireducens. The 
discovery of genes that are specifically required for the metabolism of benzene, but not 
phenol in G. metallireducens is an important step in potentially identifying the mechanisms 
for anaerobic benzene activation. 
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INTRODUCTION 

Microbial degradation of the hydrocarbons in petroleum reser- 
voirs is of interest because of its potential impact on the hydro- 
carbon composition and quality of deposits as well as reservoir 
souring (Head et al., 2010). Furthermore, anaerobic oxidation 
of hydrocarbons coupled to the reduction of Fe(III) minerals 
can have a significant influence on subsurface biogeochemistry, 
including the production of magnetite which can provide mag- 
netic anomaly signals that can aid in localizing deposits (Lovley, 
1991). 

The availability of pure cultures capable of anaerobically 
degrading alkanes and substituted aromatic hydrocarbons has 
yielded a substantial understanding of the pathways for the degra- 
dation of these constituents of crude oil (Widdel and Rabus, 
2001; Foght, 2008; Meckenstock and Mouttaki, 2011; Heider 
and Schiihle, 2013). Due to a paucity of model pure cultures, 
less is known about the anaerobic degradation of benzene, a 
significant crude oil constituent, not only because of its value 
as a petrochemical/fuel component, but also because of its 



human toxicity (Lovley, 2000; Coates et al., 2002; Vogt et al., 
2011). 

Until recently, the study of anaerobic benzene degradation 
has been limited to investigations with mixed cultures with a 
focus on the ability of anaerobes to remove benzene from con- 
taminated groundwater (Widdel and Rabus, 2001; Foght, 2008; 
Lovley et al., 2011; Meckenstock and Mouttaki, 2011; Heider and 
Schiihle, 2013). Anaerobic benzene degradation has been docu- 
mented under methanogenic conditions (Grbic-Galic and Vogel, 
1987; Weiner and Lovley, 1998; Sakai et al, 2009; Masumoto 
et al., 2012) as well as with either Fe(III) (Lovley et al., 1994, 
1996), sulfate (Lovley et al., 1995; Anderson and Lovley, 2000), 
nitrate (van der Zaan et al., 2012), or an electrode (Zhang et al., 
2010) serving as the electron acceptor. In some of these mixed 
culture studies there was evidence that the first step in anaer- 
obic benzene activation was conversion to phenol (Vogel and 
Grbicgalic, 1986; Grbic-Galic and Vogel, 1987; Weiner and Lovley, 
1998; Caldwell and Suflita, 2000), whereas in other instances 
there was evidence that benzene was first metabolized to benzoate 
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(Chaudhuri and Wiesmann, 1995; Caldwell and Suflita, 2000; 
Kunapuli et al., 2008; Abu Laban et al., 2010) or toluene (Ulrich 
et al, 2005). However, technical difficulties in working with mixed 
cultures have prevented definitive studies on the mechanisms for 
anaerobic benzene activation. 

Dechloromonas aromatica grew with benzene as the sole elec- 
tron donor in an anaerobic medium with nitrate as the electron 
acceptor (Coates et al., 2001; Chakraborty and Coates, 2005; 
Chakraborty et al., 2005). Phenol was proposed to be the first 
product of benzene activation. However, the finding that the 
oxygen in the phenol produced was not derived from water 
(Chakraborty and Coates, 2005), as well as the presence of genes 
for oxygen-dependent benzene metabolism coupled with a lack 
of genes for anaerobic metabolism of phenol or other poten- 
tial aromatic intermediates (Salinero et al., 2009), suggested that 
molecular oxygen was involved in benzene metabolism, even 
though the medium was anaerobic. It has been suggested that 
D. aromatica utilizes molecular oxygen generated intracellularly 
from nitrate for benzene activation (Salinero et al., 2009; Weelink 
et al, 2010; Meckenstock and Mouttaki, 2011; Vogt et al, 2011), 
but this possibility has yet to be experimentally verified. 

The hyperthermophilic archeon, Ferroglobus placidus was the 
first organism in pure culture unequivocally found to be capa- 
ble of anaerobically oxidizing benzene (Holmes et al., 2011). 
F. placidus oxidizes benzene to carbon dioxide with Fe(III) 
as the sole electron acceptor (Holmes et al., 2011). Benzoate, 
but not phenol or toluene, transiently accumulated during 
benzene metabolism and [ 14 C] -benzoate was produced from 
[ 14 C] -benzene (Holmes et al., 2011). Genome-scale transcrip- 
tional analysis demonstrated that during growth on benzene, 
there was an increase in transcript abundance for genes specif- 
ically involved in the metabolism of benzoate, but not phenol. 
Transcript abundance for a putative carboxylase gene was higher 
during growth on benzene vs. benzoate, suggesting a potential 
candidate enzyme for the carboxylation reaction (Holmes et al., 
2011). The lack of a genetic system for F. placidus, and the tech- 
nical difficulties of working with a hyperthermophile, has lim- 
ited further investigation of the benzene degradation pathway in 
F. placidus. 

However, it has subsequently been found that Geobacter met- 
allireducens, which can be genetically manipulated (Oberender 
et al., 2012; Tremblay et al., 2012; Smith et al, 2013), is also 
capable of anaerobically oxidizing benzene with the reduction of 
Fe(III) (Zhang et al., 2012, 2013). As previously reviewed (Lovley 
et al., 201 1), Geobacter species are thought to be important agents 
for the removal of benzene and other aromatic hydrocarbons 
from a diversity of contaminated subsurface environments in 
which Fe(III) minerals are available. Multiple lines of evidence 
demonstrated that G. metallireducens metabolized benzene via a 
phenol intermediate rather than benzoate (Zhang et al, 2013). 
For example, small amounts of phenol were detected during 
growth on benzene and ls O-labeling studies demonstrated that 
oxygen was derived from water to generate phenol. Transcripts 
for genes specifically involved in the metabolism of phenol were 
more abundant during growth on benzene than during growth 
on alternative aromatic substrates, and the deletion of the genes 
for subunits of two enzymes involved in phenol degradation 



prevented the metabolism of benzene whereas deleting genes 
specific for benzoate or toluene metabolism had no impact on 
benzene metabolism. 

The conversion of benzene to phenol is an exergonic reaction 
at pH 7 under standard conditions: 

Benzene + H 2 0 -> phenol + 2H+AG 0 ' = -14.7kJ/mol 
(calculated with data from Thauer et al., 1977) 

The exergonic hydroxylation reaction of benzene by G. metallire- 
ducens is consistent with other characterized biological anaerobic 
hydroxylation reactions such as the conversion of ethylbenzene to 
(S)-l-phenylethanol and the hydration of acetylene to acetalde- 
hyde, which are exergonic reactions performed by enzymes 
belonging to oxidoreductase families (Rosner and Schink, 1995; 
Johnson et al., 2001; Kniemeyer and Heider, 2001). 

Here we report on studies designed to identify genes encoding 
enzymes involved in the initial activation of benzene to phe- 
nol. Target genes identified from previously published (Zhang 
et al., 2013) whole-genome transcriptomic analysis were deleted 
to ascertain genes required for the metabolism of benzene, but 
not phenol, that might encode candidate enzymes for benzene 
activation. 

MATERIALS AND METHODS 
ORGANISMS AND CULTURE CONDITIONS 

The bacterial strains and plasmids used in this study are listed in 
Table SI. Geobacter metallireducens (ATCC 53774 and DSM 7210) 
(Lovley et al., 1993) was routinely cultured under strict anaerobic 
conditions at 30°C with acetate (10 mM) as the electron donor 
and Fe(III) citrate (50 mM) as the electron acceptor, as previously 
described (Lovley and Phillips, 1988). 

Metabolism of labeled compounds was investigated with 
cell suspensions as previously described (Zhang et al., 2013). 
Additions for [ 14 C] -labeling studies were 2.59 x 10 5 Bq of 14 C 
[39 [iM [UL- 14 C]benzene (2.78 x 10 9 Bq mmor 1 ; Moravek 
Biochemicals, Brea, CA, USA) or 1.85 x 10 5 Bq [U- 14 C]phenol 
(2.96 x 10 9 Bq mmor 1 ; ARC, St-Louis, MO, USA)]. 

ANALYSIS OF GENE EXPRESSION 

Whole-genome microarray analysis of gene transcript abundance 
in cells metabolizing different aromatic substrates or acetate 
have previously been described (Zhang et al., 2013; NCBI GEO 
under accession number GSE33794). The microarray results were 
analyzed with Array 4 Star (DNASTAR, Madison, WI, USA). 

MUTANT CONSTRUCTION 

Mutants were constructed as described previously (Zhang 
et al., 2013). Genomic DNA was extracted with the Epicentre 
MasterPure DNA Purification Kit (Epicentre Biotechnologies, 
Madison, WI, USA). Plasmids were extracted with the QIAprep 
Spin Miniprep Kit (Qiagen, Valencia, CA, USA). PCR amplifica- 
tion was done with the Taq polymerase (Qiagen). DNA gel purifi- 
cation was done with the QIAquick gel extraction kit (Qiagen). 
Mutant alleles were constructed by replacing the coding sequences 
with a spectinomycin resistance. Briefly, upstream and down- 
stream flanking regions of the genes to be deleted (ca. 500 bp) 
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were amplified. PCR products were mixed, digested with Avrll 
(NEB, Beverly, MA) and ligated with the T4 DNA ligase (NEB). 
The resulting construct (ca. 1 kb) was cloned into pCR2.1-TOPO 
with a TOPO TA Cloning Kit (Invitrogen, Carlsbad, CA, USA). 
The spectinomycin resistance cassette with Avrll sites at both 
ends was amplified with pRG5 (Kim et al, 2005) as a template. 
The spectinomycin resistance cassette was Avr/7-digested and lig- 
ated into the Avrll site located between the flanking regions of 
the genes to be deleted. Plasmids bearing the mutant alleles were 
linearized by restriction enzyme digestion. The linearized plas- 
mids were electroporated into G. metallireducens as described 
previously (Tremblay et al, 2012). Genotypes of the mutant 
strains were confirmed with PCR and the absence of undesired 
mutations was confirmed with Sanger sequencing. 

RESULTS AND DISCUSSION 

In order to identify genes that might be specifically associated 
with the initial conversion of benzene to phenol, gene transcript 
abundance was compared in cells metabolizing benzene and cells 
metabolizing phenol. There were 1 1 genes with higher transcript 
abundance in cells metabolizing benzene than in cells metabo- 
lizing phenol (Table 1). Each of these genes was also more highly 
expressed in cells metabolizing benzene than in cells metabolizing 
benzoate, toluene, or acetate (Table 1). 

Five of these genes with higher transcript abundance in 
benzene-metabolizing cells were predicted to be involved either 
in flagella synthesis, in the CRISPR system, or to be a toxin or a 
porin (Table 1), suggesting that they are unlikely to be involved 
in an enzymatic attack on benzene. Therefore, the function of 
these genes was not further investigated. Strains in which one of 
the remaining six genes was deleted were constructed to deter- 
mine whether the loss of the gene specifically impacted benzene 
metabolism. 

Strains lacking either Gmet 0244, Gmet 2410, or Gmet 3229, 
oxidized benzene at rates just slightly lower or comparable 
to wild type (Figure 1), suggesting that they did not encode 
enzymes important in benzene metabolism. Gmet 0244 encodes 



a hypothetical protein, whereas Gmet 2410 and Gmet 3229 are 
predicted to encode a superoxide dismutase, and a thioredoxin- 
related protein, respectively. The strain in which Gmet 2833 was 
deleted was deficient in phenol metabolism as well as benzene 
metabolism (Figure 2), indicating that this gene of unknown 
function did not have a role specific to benzene metabolism. 

The strain lacking a functional Gmet 3376 oxidized benzene 
slightly slower and to a lesser extent than wild-type, whereas this 
strain metabolized phenol just as well as wild-type (Figure 2). 
Gmet 3376 encodes a hypothetical protein of 219 amino acids 
with no known homolog in any other Geobacter species. The abil- 
ity of the Gmet 3376-deficient strain to continue to metabolize 
benzene at substantial rate suggests that it does not encode an 
enzyme that is required for benzene activation. 

In contrast, deletion of the monocistronic gene Gmet 0232 
specifically inhibited the metabolism of benzene, but not phenol 
(Figure 2). Gmet 0232 encodes for a hypothetical protein of 281 
amino acids with homologs found only in Geobacteraceae family 
members and in Nitrosomonas sp. (strain Is79A3). 

The transcript abundance of the adjacent gene, Gmet 0231 
(Figure 3), was elevated in benzene-metabolizing cells (44.6- 
fold compared to acetate-metabolizing cells; P-value = 0.00004), 
but the difference in transcript abundance between benzene- 
and phenol-metabolizing cells was too low to be considered to 
be significantly higher in benzene-metabolizing cells (1.7-fold 
higher; P-value = 0.003). Gmet 0231 is predicted to be a mono- 
cistronic gene encoding a zinc-containing oxidoreductase of the 
NADPH:quinone oxidoreductase family. A strain in which Gmet 
0231 was deleted had a phenotype similar to the strain deficient 
in Gmet 0232; phenol was metabolized as well as in wild-type, but 
benzene metabolism was inhibited (Figure 2). 

IMPLICATIONS 

The results demonstrate that G. metallireducens specifically 
requires genes Gmet 0231 and 0232 for the metabolism of ben- 
zene, but not phenol. This raises the possibility that Gmet 023 1 or 
Gmet 0232, or both, encode an enzyme(s) involved in the initial 



Table 1 | Genes up-regulated at least two fold in G. metallireducens when benzene is the electron donor (P-value cutoff < 0.1; NCBI GEO under 
accession number GSE33794). 
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Annotation 
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Fold change: benzene vs. 










Phenol 


Benzoate 


Toluene 


Acetate 


Gmet 0232 a 


Conserved hypothetical protein 




3 


20.2 


2.6 


18.1 


Gmet 0244 a 


Conserved hypothetical protein 




3.9 


14 


5 


19.5 


Gmet 0802 


Putative porin 




2.4 


9.9 


3.5 


14.9 


Gmet 1065 


CRISPR-associated protein 


casD 


2.3 


13.3 


2.6 


13 


Gmet 2410 a 


Superoxide dismutase 


sodA 


6.1 


2.3 


2.7 


4.3 


Gmet 2833 a 


Conserved hypothetical protein 




2.7 


5.5 


3 


5.1 


Gmet 3104 


Flagellar operon protein of unknown function DUF3766 




3.2 


8.4 


9 


8 


Gmet 3105 


flagellar hook capping protein 


flgD 


2.5 


5.9 


7.4 


5.4 


Gmet 3076 


Toxin, PIN family 




2.7 


12.4 


3.4 


9.5 


Gmet 3229 a 


Thioredoxin-related protein disulfide reductase, putative 




2.1 


6.5 


2.9 


4.4 


Gmet 3376 a 


Conserved hypothetical protein 




2.2 


137.3 


3.6 


141.6 
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activation of benzene to produce phenol, or serve an important 
ancillary role in this phase of benzene metabolism. 

Gmet 0232 homologs are found in two Geobacter species, 
Geobacter sulfurruducens (Caccavo et al., 1994), and Geobacter 
uraniireducens (Shelobolina et al., 2008), that are unable to 
anaerobically oxidize monoaromatic hydrocarbons. Gmet 0231 
homologs are found in multiple species throughout the three 
kingdoms of life that are not known to degrade monoaromatic 
hydrocarbons, including several Geobacter species. However, the 



A 

100 



B 




Days 

FIGURE 1 | Production of 14 C0 2 by cell suspension of strains lacking 
Gmet 0244, Gmet 2410, or Gmet 3229 from (A) [ 14 C]benzene or (B) 
[ 14 C]phenol. The results are the mean and standard deviation for triplicate 
cell suspensions. 



co-localization of Gmet 0231 and Gmet 0232 found in G. met- 
allireducens has not been observed in any other microbial 
genome sequence, suggesting that a functional relation between 
the proteins encoded by Gmet 0231 and Gmet 0232 might 
be necessary for benzene conversion to phenol by G. metal- 
lireducens. Now that Gmet 0231 and Gmet 0232 have been 
identified as having an important role in the initial steps in ben- 
zene metabolism, biochemical studies of the gene products is 
warranted. 




Zinc-containing 
oxidoreductase 
DoxXfamily; related to 
putative NADPHiquinone 
oxidoreductaseoxidoreductase 



Hypothetical 
protein 



Hypothetical 
proteins 



Hypothetical 
proteins 



^0230 0231 0232 0233^023 4^ 0235 ^ 023 6^> 



FIGURE 3 | The Gmet 0232 genomic region. 
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Also warranted are studies on the function of genes localized 
near Gmet 0231 and Gmet 0232 in the G. metallireducens genome. 
For example, Gmet 0230 is predicted to encode an oxidoreductase 
of the DoxX family, which has sequence similarities with DoxD, a 
subunit of the membrane-bound thiosulphate:quinone oxidore- 
ductase (Muller et al., 2004). Genes adjacent to Gmet 0232 are 
predicted to encode hypothetical proteins of unknown function 
(Figure 3), suggesting that genetic analysis of their function in 
benzene metabolism may also be warranted. 

It has been suggested (Coates et al, 2002) that enzymes 
with functions similar to the ethylbenzene dehydrogenase found 
in Aromataleum aromaticum (Johnson and Spormann, 1999; 
Kniemeyer and Heider, 2001) or the acetylene hydratase of 
Pelobacter acetylinus (Schink, 1985; Rosner and Schink, 1995) 
might be involved in the anaerobic hydroxylation of benzene. 
None of the genes investigated here have any homology to the 
components of either of those enzyme complexes. Thus, these 
results suggest that G. metallireducens may possess a novel bio- 
chemistry for the anaerobic hydroxylation of benzene. 
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